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Abstract 

We construct a chiral Lagrangian for doubly heavy baryons and heavy mesons that is invariant 
under heavy quark-diquark symmetry at leading order and includes the leading ©(l/mg) symmetry 
violating operators. The theory is used to predict the electromagnetic decay width of the J = ^ 
member of the ground state doubly heavy baryon doublet. Numerical estimates are provided for 
doubly charm baryons. We also calculate chiral corrections to doubly heavy baryon masses and 
strong decay widths of low lying excited doubly heavy baryons. 
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Heavy quark-diquark symmetry relates mesons with a single heavy quark to antibaryons 
with two heavy antiquarks. Savage and Wise [l| argued that quark-diquark symmetry 
was realized in the heavy quark limit of Quantum Chromodynamics (QCD) and studied 
th,s symmetry us„>g the methods of Heavy Quark Effective Theory (HQET) g Reeently, 
Refs. j3|, 1^ derived effective Lagrangians for heavy diquarks within the framework of Non- 
Relativistic QCD (NRQCD) 0,0,01 • These papers obtain a prediction for the hyperfine 
splitting of the ground state doubly heavy baryons in terms of the ground state heavy 
meson hyperfine splitting^. Heavy quark-diquark symmetry also relates other properties of 
heavy mesons and doubly heavy baryons. A useful tool for studying low energy strong and 
electromagnetic interactions of heavy hadrons is heavy hadron chiral perturbation theory 
(HHxPT) 0, y,0. This theory has heavy hadrons, Goldstone bosons, and photons as its 
elementary degrees of freedom and incorporates the approximate chiral and heavy quark 
symmetries of QCD. In this paper we derive a version of HH^PT that includes doubly 
heavy baryons and incorporates heavy quark-diquark symmetry. The theory is used to 
calculate chiral corrections to doubly heavy baryon masses and to obtain model-independent 
predictions for the electromagnetic decay of the </ = | member of the ground state doubly 
heavy baryon doublet. Our formulae are applicable to either doubly bottom or doubly charm 
baryons, and we give numerical estimates for the case of doubly charm baryons. We also 
discuss the low lying excited doubly heavy baryons, show how these states can be included 
in the effective theory, and calculate their strong decay widths. 

Motivation for this work comes from the SELEX experiment's recent observation of states 
which have been tentatively interpreted as doubly charm baryons 12 . and also the 



COMPASS experiment, which in its second phase run in 2006 hopes to observe doubly 
charm baryons Many aspects of the SELEX states are difficult to understand. States 
observed by SELEX include the S+(3520), which decays weakly into A+vr+i^" 11] as well 
as pD^K^ [13!], the E!++(3460), which decays weakly into A'^K^n^'n-^ 12], and a broader 
state, S++(3780), also seen to decay into A+i^^vr+vr"'" [l^. The ground states of the S+ 
and S^"*" are related by isospin symmetry and therefore should differ in mass by only a 
few MeV, so the observed difference of 60 MeV seems implausible. On the other hand, an 



^ The formula for the hyperfine sphttings in Ref. jlj differs from the correct formula in Refs. [3|, ^ by a 
factor of 2. 
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unpublished talk [3] and conference proceedings |3] present evidence for additional states, 
S+(3443) and 2++ (3541). If these states exist the isospin splittings are closer to theoretical 
expectations, but still quite large. The difference between the mass of the S+(3520) and 
H+(3443) is 77 MeV, and the splitting between S++(3541) and S++(3460) is 81 MeV. These 
splittings agree remarkably well with calculations of the doubly charm hyperfine splittings 



in quenched lattice QCD 
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ISl |l9| and are within ~ 25 MeV of the heavy quark-diquark 
symmetry prediction obtained in Refs. 0,0], an acceptable discrepancy given the expected 
0(Aqcd/'^c) corrections. However, interpretation of the S+(3520) as the ^ = | member 
of the ground state doublet is impossible to reconcile with the fact that the S^(3520) 
is observed to decay weakly because if the S+(3520) is not the ground state of the ccd 
system it should decay electromagnetically. There are also discrepancies between the weak 



decay lifetimes predicted by HQET j20|, |21|, 22|(~ 100 fs) and the observed lifetimes (< 
33 fs) [ill Q]. Production cross sections are also poorly understood within perturbative 
QCD |23i. However, the SELEX states are observed in the forward region, (xf) ~ 0.3, 

iTnn 

where nonperturbative production mechanisms such as intrinsic charm j2j, |2^ |2a] or parton 
recombination 2^ may be important. 

Even if there is difficulty interpreting the SELEX data, doubly charm baryons must exist 
and are expected to have masses of approximately 3.5 GeV, where the SELEX states are. 
In light of existing and future experimental efforts to observe doubly charm baryons, it is 
desirable to have model independent predictions for other properties besides the relation 
for the hyperfine splittings derived in Refs. Therefore it is important to develop 

theoretical tools for analyzing the properties of doubly heavy baryons systematically. 

Savage and Wise |1| wrote down a version of heavy quark effective theory (HQET) which 
includes diquarks as elementary degrees of freedom and derived a formula relating heavy 
meson and doubly heavy baryon hyperfine splittings. HQET only separates the scales Aqcd 
and rriQ, where rriQ is the heavy quark mass. The dynamics of a bound state of two heavy 
quarks is characterized by additional scales rnqv and rnqv'^, where v is the typical velocity 
of the heavy quarks in the bound state. The correct effective theory for hadrons with two 
heavy quarks is Non-Relativistic QCD (NRQCD) 5], which properly accounts for the scales 
rnqv and rnqv"^. Analysis of heavy diquarks within the framework of NRQCD was recently 
performed in Refs. j^l |^. These papers derived Lagrangians for diquark fields starting 
from NRQCD and obtained the correct heavy quark symmetry prediction for the hyperfine 
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splittings of the doubly heavy baryons. For simplicity, we will consider only one flavor of 
heavy quark. The lowest mass diquark will consist of two heavy antiquarks in an orbital 
S-wave in the 3 representation of color SU{3). Then Fermi statistics demands that they 
have total spin one. In the rest frame of the heavy quark and lowest mass diquark, the 
Lagrangian to ©(l/mg) is Q] 

C = {tDo- ] h + V^ ■ (tDo + 6 - —] V 




2mQ 2 2mQ 2 ^ ^ 

Here h is the heavy quark field, V is the field for the diquarks, the A"/2 are the SU{3) 
color generators, Tr[A"A^] = 2 5"^, Dq and D are the time and spatial components of the 
gauge covariant derivative, respectively, i?" is the chromomagnetic field, and rriQ is the 
heavy quark mass. The term proportional to 6 is the residual mass of the diquark. The 
heavy antiquarks in the diquark experience an attractive force and therefore the mass of 
the diquark is not 2mQ but 2mQ — 6, where 6 is the binding energy. This residual mass 
can be removed by rephasing the diquark fields. Physically, this corresponds to measuring 
diquark energies relative to the mass of the diquark, rather than 2mQ. Once this is done the 
Lagrangian, at lowest order in I/itlq, is invariant under a U{5) symmetry which permutes 
the two spin states of the heavy quark and the three spin states of the heavy antiquark. The 
f/(5) symmetry is broken by the 0{l/mQ) kinetic energy and chromomagnetic couplings of 
the heavy quark and diquark. The latter terms are responsible for the hyperfine splittings. 

The ground state doublet of heavy mesons is usually represented in HHxPT as a 4 x 
4 matrix transforming covariantly under Lorentz transformations, and transforming as a 
doublet under SU{2) heavy quark spin symmetry, 

Hv = (^) {P7l, - . (2) 



Here P*^ is the = 1^ vector heavy meson field which obeys the constraint v^P*'^ = 0, 
where f is the four-velocity of the heavy meson. is the = 0~ pseudoscalar heavy 
meson field. The superfield obeys the constraints ■^H^ = —H^fi = H^, so only has 
four independent degrees of freedom. These can be collected in a 2 x 2 matrix. For example, 
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in the heavy meson rest frame where v^^ — (1, 0, 0, 0), 

B. J" (3) 
y 

where we have used the Bjorken-Drell conventions for 7^ and 75. For a process such as 
the weak decay B — > Dii/, in which the initial and final heavy hadrons have different 
four-velocities, the covariant representation of fields is needed to determine heavy quark 
symmetry constraints on heavy hadron form-factors. However, for studying low energy 
strong and electromagnetic interactions in which the heavy meson four-velocity is conserved 
(up to 0(Aqcd/^q) corrections), it is also possible to work in the heavy meson rest frame 
and use 2x2 matrix fields. This makes some calculations simpler and we find it easiest to 
formulate the extension of HHxPT with U (5) quark-diquark symmetry in this frame. We 
define the heavy meson field in our theory to be 

Ha^P:-a + Pa, (4) 

where a is an SU (3) flavor anti-fundamental index and the a are the Pauli matrices. Since 
we have chosen to work in the heavy meson rest frame, Lorentz covariance is lost and 
the symmetries of the theory are rotational invariance, SU (2) heavy quark spin symmetry, 
parity, time reversal and SUl{^) x SUr{3) chiral symmetry. Under these symmetries the 
field Ha transforms as 



rotations 


K - 


= UHJJ^ 


heavy quark spin 


K - 


- SHa 


parity 


K - 


= -Ha 


time reversal 


K - 


- —a2H*a2 


SUl{S) X SUr{3) 


K - 


-- H,Vl . 



(5) 

Here U and S are 2x2 rotation matrices and V^^ is an SU (3) matrix which gives the 
standard nonlinear realization of SUl{2>) x SUr{2>) chiral symmetry. In the two component 
notation the HHxPT Lagrangian is: 

C = TY[Hl{iDo)baH,] - gTx[HlH, a ■ A,a] + ^Tr[i7t Ha a'] . (6) 

The last term breaks heavy quark spin symmetry and Ajj is the hyperfine splitting of the 
heavy mesons. The time component of the covariant chiral derivative is (£'o)6a) ^60 is the 



spatial part of the axial vector field, and g is the heavy meson axial coupling. Our definitions 
for the chiral covariant derivative, the axial current, and the Lagrangian for the Goldstone 
boson fields are the same as Ref. j29[. 

We are now ready to generalize the Lagrangian to incorporate the doubly heavy baryons 
and the U{5) quark-diquark symmetry. The field Ha transforms like the tensor product of 
a heavy quark spinor and a light antiquark spinor. (This is how representations of heavy 



hadron fields were constructed in Ref. [30|.) Writing the field with explicit indices, {Ha)a(3, 
the index a corresponds to the spinor index of the heavy quark and the index (3 is that of 
the light antiquark spinor. In the theory with quark-diquark symmetry, the heavy quark 
spinor is replaced with a five-component field, the first two components corresponding to 
the two heavy quark spin states and the last three components corresponding to the three 
spin states of the diquark: 




= • (7) 



In terms of the kinetic terms of the Lagrangian in Eq.(^ are 

C = QltDoQ^. (8) 

The fields in HHxPT with heavy quark-diquark symmetry transform as tensor products of 
the five component field and a two-component light antiquark spinor. Thus, the 2x2 
matrix field Ha is promoted to a 5 x 2 matrix field 

Ha,al3 — ^ '^a,^J.|3 = Ha^afi + Ta^ip . (9) 

Here the index /i takes on values between 1 and 5, a,/? = 1 or 2, and i = 3,4, or 5. The 
doubly heavy baryon fields are contained in Ta^j/j. Under the symmetries of the theory Tia 
transforms as 

rotations H'a = TZHaU^ 
heavy quark spin = SHa 

parity H'a = -Ha 
time reversal H'a = —T,2Hla2 
SUj^{3) X SUn{3) K = H,Vl. (10) 
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The matrix S is now an element of U{5) and 7^ is a 5 x 5 reducible rotation matrix 

n,.-("-°], (u, 

where Uaf3 is an SU (2) rotation matrix and Rij is an orthogonal 3x3 rotation matrix related 
to U by U'^aiU = R-ijUj- The 5x5 matrix appearing in the time reversal transformation is 



= . (12) 




Under rotations the field Ta^ip transforms as T^^^^ = RijTaj-yU^^p. Ta^ip can be further 
decomposed into its spin-| and spin-^ components, 



Ta,il3 — (^*a,ip + ^'^'''T ^iP^ 



(13) 



where S* .^^ and S^^^ are the spin-| and spin-| fields, respectively. The factor of \/2 is 
a convention that ensures that the kinetic terms of the doubly heavy baryon fields have 
the same normalization as the heavy meson fields. The field S*^^ obeys the constraint 

The f/(5) invariant generalizations of the first two terms of Eq. (jH)) are simply obtained 
by making the replacement Ha 'Ha- To determine the proper generalization of the f/(5) 
breaking term we note that the chromomagnetic couplings in Eq. can be written as 

^fiiS„..B.fs„, (14) 

where the S^j, are the 5x5 matrices 



= , (15) 




and {T^)jk = —i^ijk- It is now obvious that the correct generalization of Eq. (jHI) is 

£ = TilHliiDoUnb] - gTilHlHba ■ A,a] + ^TiiHlE^Haa^ 
= TT[Hl{tDo)baHb] - g Tt[HIH, a ■ A,a] + ^Tt[HI a' Ha a'] 

+TT[T^{zDo)baTb] - g TtMT, a ■ A,a] + ^Tr[Tt T' a'] . (16) 
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The last line of Eq. ()16p contains the terms relevant for doubly heavy baryons. Heavy 
quark-diquark sjTiimetry relates the couplings in the doubly heavy baryon sector to the 
heavy meson sector. The propagator for the spin-| doubly heavy baryon is 

2{ko + AH/2 + ie) ' 
while the propagator for the spin-| doubly heavy baryon is 



2{ko- An/A + te) 2{ko - Ah + ie) 

The projection operator Via,j/3 satisfies cr^„ 'Pia,j/3 = 'Piaj/s c"^^ = 0. Comparison of the poles 
of the propagators shows that the hyperfine splitting for the doubly heavy baryons is ^A^, 
reproducing the heavy quark-diquark symmetry prediction 

3 

rris* — m-s = -{nip* — mp) , (17) 

obtained in Refs. 0,0]. 

We can also include operators which mediate electromagnetic decays. The Lagrangian 
for electromagnetic decays of the heavy mesons in the two-component notation is 31 1 

C = ^TT[HlHb a ■ B Q,b] + -^Q'Tr [i/tc? . BH^] , (18) 
2 2mQ 

where Qab = diag(2/3, — 1/3, — 1/3) is the light quark charge matrix, P is the parameter 
introduced in Ref. j^ll, and Q' is the heavy quark charge. For charm, Q' = 2/3. The first 
term is the magnetic moment coupling of the light degrees of freedom and the second term 
is the magnetic moment coupling of the heavy quark. Both terms are needed to understand 
the observed electromagnetic branching fractions of the D*~^ and Z)*° 0]. The magnetic 
couplings of the heavy quark and diquark are 

= g'/it ff.Bh Q' -Bx V 

2mQ niQ 

■Q'Qlt'^^-BQ,, (19) 



2mr 



where the S^^, are the 5x5 matrices 



\ |. (20) 



-2% 



jk 



The magnetic coupling of the diquark has the opposite sign as that of the heavy quark 
because it is composed of two heavy antiquarks. The coefficient of the chromomagnetic 
couphng of the diquark in Eqs. (jlll4p is a factor of 2 smaller than the coefficient of the 
electromagnetic coupling of the diquark in Eq. (jl9|) due to a color factor. The magnetic 
couplings in the HH^PT Lagrangian for heavy mesons and doubly heavy baryons are 

L = ^Trinl n,a-B Q^b] + -^Q'Tr[7^1 S' -BU,]. (21) 
Z ZniQ 

The part of this Lagrangian involving the doubly heavy baryon fields is 

C = ^Tr [Ti T,a-B Q„,] - — Q'Tr [T] f-BT,]. (22) 

2 niQ 

This can be used to obtain the following tree level predictions for the electromagnetic decay 
widths: 

3 V t^qJ "^p* 

W: - 2.7] = ^ (pQaa - ^) ' ^i?? . (23) 

Here is the photon energy. These results could also be obtained in the quark model, 
with the parameter j3 = l/mg, where rriq is the light constituent quark mass. The effective 
theory allows one to improve upon this approximation by including corrections from loops 
with light Goldstone bosons, which give 0{y/fn^) corrections to the decay rates j3l|. If these 
loop corrections are evaluated in an approximation where heavy hadron mass differences are 
neglected, the correction to the above formulae can be incorporated by making the following 
replacements j^^] 

For charm mesons, hyperfine splittings are ~ 140 MeV and the SU{3) splitting is ~ 100 
MeV, while for bottom mesons the hyperfine splittings are ~ 45 MeV and SU{3) splitting 
is ~ 90 MeV. The approximation of neglecting heavy hadron mass differences and keeping 
Goldstone boson masses is reasonable for kaon loops but not for loops with pions. However, 



the largest 0{y/m^) corrections come from loops with kaons. When data on double heavy 
baryon electromagnetic decays is available, more accurate calculations along the lines of 
Ref. should be performed. In this paper, we will use Eqs. and to obtain 
estimates of doubly charm baryon electromagnetic decay widths. 

Currently r[D*"^] is measured to be 96 ± 22 keV, while there is only an upper limit for 
T[D*^]. The branching ratios for the D*+ decays are Bt[D*+ D^7t+] = 67.7 ± 0.5%, 
Br[D*+ D+ttO] = 30.7 ± 0.5% and By[D*+ D+7] = 1.6 ± 0.4%. The branching ratios 
for D*^ decays are Bt[D*^ D^n^] = 61.9 ± 2.9% and Bi[D*^ D^-f] = 38.1 ± 2.9%. 
Isospin symmetry can be used to relate the strong partial width of the D*^ to the known 
strong partial width of the D*^. Then the measured branching fractions of the D*^ can be 
used to obtain the partial electromagnetic width of the D*^. We find 

r[L)*o^D%] = 26.1±6.0keV 

r[D*+^D+7] = 1.54 ± 0.35 keV, (25) 

where the error is dominated by the uncertainty in r[D*+]. r[D*^ D^'-f] is suppressed 
because of a partial cancellation between the magnetic moments of the light degrees of 
freedom and the charm quark. Using the partial widths in Eq. and the formulae in 
Eqs. fl23j) and (j2H), we obtain predictions for doubly charm baryon electromagnetic decays 
in Table |l] 

In our calculations of the doubly charm baryon decay widths the factor ms/ms* in 
Eq. ()23|) has been set equal to one. For the expected masses and hyperfine splittings of 
the doubly charm baryons, this factor changes the predictions for the widths by less than 
3%. The fits are labeled in the left hand column of Table d In the fits labeled QM we 
have not included the O(y^m^) corrections in Eq. ()24j) . Therefore, these predictions for the 
doubly charm baryon electromagnetic decays are the same as what would be obtained in 
the quark model. The values of the parameters f3 and for each fit are shown along with 
the predictions for the electromagnetic decay widths. In QM 1, we have treated (3 and rric 
as free parameters and fit these to the central values in Eq. (j25|) . In QM 2 we have set 
TTLc = 1500 MeV and performed a least squared fit to (3. In the fits labeled xPT, we have 
included the leading 0{y/m^) chiral corrections in Eq. (j21I). We have used = 130 MeV, 
Jk = 159 MeV, and g = 0.6 which is extracted from a tree level fit to the D*^ width. In 
xPT 1, we fixed /3 and rric to reproduce the central values in Eq. (|^. In xPT 2, we set 
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Fit 



mc(MeV) 



r[H:++] (keV) 



mt] (keV) 



QM 1 
QM 2 
XPT 1 
XPT2 



379 



356 



272 



276 



1863 



1500 



1432 



1500 



3.3 



80 MeV 



2.6 ' 



3.4 



80 MeV 



3.2 



80 MeV 



80 MeV J 



^•^ ' 80 MeV ) 



^■^ ' 80Mev) 



2.3 



80 MeV 



3.3 



V 

80 MeV y 



TABLE I: Predictions for the electromagnetic widths of the S*^ and 
in the text. 



The fits are explained 



rric = 1500 MeV and performed a least squares fit to f3. There are several sources of error in 
the calculation. We expect 30% theoretical errors due to heavy quark symmetry breaking 
effects, 30% errors due to higher order SU{3) breaking effects, and 25% uncertainty from 
the experimentally measured value of r[D*+] leading to at least 50% error in the predictions 
in Table EI 

Chiral perturbation theory and the nonrelativistic quark model give similar size estimates 
for the S*^ electromagnetic decay widths which are expected to be ~2-3 keV if the hyperfine 
splitting is 80 MeV. The electromagnetic decay should completely dominate any possible 
weak decay, even if the weak decay rates are an order of magnitude greater than calculated 



in Refs. |20l . 
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22j . The quark model predicts r[S*+"'"] slightly greater than r[S*+]. This 



is in contrast with the charm meson sector where the magnetic moment of the light degrees 
of freedom and the magnetic moment of the charm quark add constructively to give a large 
r[D*^ — s> -D°7] and destructively to give a small r[D*^ D^l]- In the doubly heavy 
baryon sector, the relative sign of the magnetic moments is reversed, and both decay rates 
are approximately the same. In fact from Eq. (j^Hj) . we can see that for (3 = A/rric the 
two rates are exactly equal in the quark model. Fits to the D* electromagnetic decays 
yield values of P and nic that are close to this point in parameter space. Including the 
0{y/ffhi) corrections from chiral perturbation theory, the most important effect is the kaon 
loop correction whose contribution to the decay has opposite sign as the contribution 
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from (5 at tree level, therefore suppressing the S*++ decay relative to S*+. 

The theory can also be used to compute chiral corrections to doubly heavy baryon masses. 
The one loop corrections to the hadron masses are 

.2 /5 



'—'a 



i,b 



9 



167r2/f V9 



) 4 \ 



167r2/? V9 



1 8 



9 



y^.^ab 



9 



Lb 



EC, 



Lb 



"''167r2/2 



K{mH* -mH,,mi,fi) 



1 2 

-K{mH, - mH*,mi, n) + -K{mH* - mH*,mi, fi) 



(26) 



Here rrii and /j are the mass and decay constant of the Goldstonc boson in the one loop 
diagram and is a factor which comes from SU{2>) Clebsch-Gordan coefficients in the 
couphngs. For loops with charged pions we have — Cjf = 1, for loops with neutral 
pions Cfi = = \i for loops with kaons C|[ = = 1 (i =1 or 2), and for loops with 77 
mesons C^^ — C22 — \ and C33 = |. The function K{5,m,iJi) is related to the finite part of 
the integral 



1 



1 



{2T:)^l'^-ml + ielQ-5 + ie {An) 



(27) 



evaluated using dimensional regularization in the MS scheme. We find 

.2\ 



K{5, m, = (-2 5^ + 3 5) In 
where 



m 



+ 25 {5^- m') F 



m 



+ A5-'-h5m\ (28) 



F{x) = 2 



TT 



- Tan-1 



= -2 



X 



X 



x^ 



X 



In ^x + Vx"^ — 1 ^ 



< 1 

|x| > 1 , 



and fi is the rcnormalization scale. The fi dependence in the one loop calculation is cancelled 
by counterterms that have not been included. 

We are interested in how the one loop corrections affect the leading order prediction for 
the hyperfine splittings. Unfortunately, it is impossible to give a reliable estimate without 
knowing the numerical value of the counterterms required to cancel the /i dependence in 
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the nonanalytic contribution. Furthermore, to compute the contribution from kaon loops, 
one must know the masses of doubly charm strange baryons which have not been observed. 
We will assume that the ground state doubly charm strange baryons are 100 MeV higher in 
mass than their nonstrange counterparts, similar to the D meson system. We work in the 
isospin limit and use g = 0.6, = 140 MeV, = 137 MeV, mx = 496 MeV, = 548 
MeV and the experimental values of the pseudoscalar meson decay constants: /■„ = 130 
MeV, fx = 159 MeV, and = 156 MeV. The nonanalytic part of the one loop correction 
to the nonstrange doubly charm baryon hyperfine splitting is 

-7.0 MeV /i = 500 MeV 
Smp,^-6m^_=l 8.1 MeV /i = 1000 MeV , (29) 

16.9 MeV /i = 1500 MeV 

where we have shown our results for three values of fi. For these choices of /i the nonanalytic 
part of the chiral correction varies between -7 MeV and +17 MeV. The nonanalytic part of 
the chiral correction to the doubly charm baryon hyperfine splitting is quite sensitive to the 
choice of fi, and lies within 15% of the tree level prediction. We also calculate the correction 
to the hyperfine splitting relationship of Eq. ()17|) and find for the masses in the nonstrange 
sector 
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3.9 MeV /i = 500 MeV 

5.3 MeV /i = 1000 MeV • (30) 

6.1 MeV /i = 1500 MeV 



The nonanalytic correction to the symmetry prediction is small (< 10 MeV) and relatively 
insensitive to the choice of /i. Chiral perturbation theory predicts the nonanalytic depen- 
dence of the doubly heavy baryon masses on the light quark masses, and generalized to 
include the effects of quenching as well as other lattice artifacts, formulae such as those 
in Eq. should be useful for chiral extrapolations of doubly heavy baryon masses and 
hyperfine splittings in lattice simulations. 

Finally, we discuss excited doubly heavy baryons. There are two types of excitations in 
the doubly heavy baryon system: excitations of the light degrees of freedom and excitations 
of the diquark. Excitations of the first type are related to analogous excitations in the heavy 
meson sector by heavy quark-diquark symmetry. The lowest lying excited charm mesons 
are in a doublet of = 0+ and 1^ mesons with masses approximately 425 MeV above 
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the ground state in the nonstrange sector 



32 



mm 

and 350 MeV above the ground state 



in the strange sector [35|,|36[. In the nonstrange sector these states decay via S-wave pion 
emission and have widths in the range 250-350 MeV, while in the strange sector the strong 
decay is via tt^ emission which violates isospin, and therefore the states are very narrow 
with widths expected to be of order 10 keV j37|. These states have light degrees of freedom 
with angular momentum and parity = The doubly charm baryons related to the 
even-parity excited charm mesons by quark-diquark symmetry are a doublet with = |^ 
and = 1^. The excited charm mesons and doubly charm baryons can be incorporated 
into HH^PT with a 5 x 2 matrix field which is like the field Tif^jj except 5^/3 has 
opposite parity. The excitation energies and strong decay widths of these excited doubly 
charm baryons should be similar to their counterparts in the charm meson sector. Since 
the excited S++(3780) state observed by SELEX is only 320 MeV above the S++(3460), the 
lowest mass S'^^ candidate, and its width is considerably less than 300 MeV, it does not 
seem likely that this excited doubly charm baryon is related to the excited charm mesons 
by heavy quark-diquark symmetry. 

This is not unexpected as the lowest lying excited doubly charm baryons are not excita- 
tions of the light degrees of freedom but rather states in which the diquark is excited. The 
lowest mass excited diquark is a P-wave excitation. Because of Fermi statistics the diquark 
is a heavy quark spin singlet. The diquark's orbital angular momentum couples with the 
angular momentum of the light degrees of freedom to form baryons with = |^ and 
7^ = 1'*', which we will refer to as and S^*, respectively. The next lowest lying states 
are doubly heavy baryons with a radially excited diquark, which form a heavy quark doublet 
with = \ and = ^ baryons, which we will refer to as and S^*, respectively. If 
the heavy antiquarks are sufficiently heavy that the force between them is approximately 
Coulombic, they interact via a potential which is 1/2 as strong as the potential between the 
quark and antiquark in a quarkonium bound state. Therefore we expect the excitation ener- 
gies of the charm diquarks to be significantly smaller than the analogous excitation energies 
in charmonium. Quark model calculations of excited doubly charm baryons predict that the 
and H^* states are about 225 MeV above the and S*^, respectively, and that that 
the heavy quark doublet containing and S^* is about 300 MeV above the ground state 



doublet 



38, 



3S 



40 



41 



42]. These excitation energies are about 1/2 the corresponding exci- 



tation energies in the charmonium system: rrih^ — mj/^ = 430 MeV and m^i — mj/^ = 590 



14 



MeV. The charm diquark excitation energies are less than the expected excitation energy 
of the hght degrees of freedom and therefore the lowest lying excited doubly charm baryons 
have excited diquarks. Excitation energies of a diquark made from two bottom quarks are 
similar to the excitation energies of a diquark made from charm, so the same conclusion 
holds for doubly bottom baryons. 

The doubly heavy baryons with P-wave excited diquarks decay to the ground state via 
S-wave pion emission. These decays violate heavy quark spin symmetry because the total 
spin of the diquark is changed in the transition. The Lagrangian for the excited and S^* 
states, including kinetic terms, residual mass terms and terms which mediate the S-wave 
decays, is 

+2 Xy, (Ht EUl + h.c. ) + 2 A3/2 El*Al + h.c. ) . (31) 

The strong decay widths of the P-wave excited nonstrange doubly charm baryons are 

/I \ 



/I 

T[El ^ S,, ^] = ^ l^-^KoP.o + El^p^^^ = A?/2 111 MeV . 
To obtain numerical estimates, we have assumed the masses m= = 3440 MeV, m=. = 3520 

' ^CC I ^cc 

MeV, Tn^^v^ = 3665 MeV and = 3745 MeV, corresponding to a diquark excitation 

energy of 225 MeV. We sum over both charged and neutral pion decay modes. The coupling 
constants A1/2 and A3/2 are 0(Aqcd/'^q) so we should expect this suppression makes A1/2 
and A3/2 < 1. Therefore these states could be quite narrow despite decaying via S-wave pion 
emission. The small widths are due to the small excitation energy which leaves little phase 
space for the decay. If the excitation energy is increased to 280 MeV, the widths are twice 
as large. Like the isospin violating decays of the D* [4^ and the even-parity excited Dg 
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4j|, the excited doubly heavy strange baryons below the kaon threshold decay 



through a virtual rj which mixes into a vr". Denoting the ground state doubly charm strange 
baryons as fii*^ and the P-wave excited doubly charm strange baryons as Q.^c*^ we obtain 
the following formulae for the isospin violating strong decay widths 

15 



Here 9 = 0.01 is the tt^ — r] mixing angle. We expect these widths to be in the range 1-5 
keV, but without knowing the masses of the Vtic and Vtcc states or the couphngs A1/2 and 
A3/2 we cannot make more precise predictions. 

The = \ and = \ doubly heavy baryons with radially excited diquarks are 
members of a heavy quark doublet we will denote whose definition in terms of component 
fields is identical to Eq. (fT^ . The Lagrangian describing this field, including terms which 
mediate its decay to the ground state, is 

C = TT[Z\{zDoU-ST'6^,)n]-gTT[Z^na-A,,] + ^TT[T^^ 

-~9 (Tr[riT,'a-l,J+h.c.) . (32) 

In the limit of infinite heavy quark mass, the light degrees of freedom in the radially excited 
doubly heavy baryons are in the same configuration as the ground state. Therefore, they 
are also related to the heavy meson ground state doublet by heavy quark-diquark symmetry. 
The axial coupling and hyperfine splitting of are the same as T^, as long as the spatial 
extent of the excited diquark, which is of order l/(mQf), is much smaller than 1/Aqcd- 
This is valid in the heavy quark limit, but could receive significant corrections in the charm 
sector. The last term in Eq. (jH^ mediates P-wave decays from the excited = | and 
= i doubly heavy baryons to the ground state. The partial decay widths are 



r[-. - ^ — -b.i r[^„ ^^,n] = c., - ^ — -b.i 

rr-' 1 r ^ 9^ 1 |3 r\-' - ^ r ^ 9^ \ \3 

r[-.--,vr] =C„,-^— n-.--.-]=C..-^^H . (33) 

Here Cab is an SU (3) factor which is 1/2 for decays involving vr" and one for decays involving 
charged pions. The radially excited doubly heavy strange baryons should also be below 
the threshold for decays into kaons, and therefore should be quite narrow. The formulae 
in Eq. ()33p can be used to obtain these decay widths as well. The isospin violating strong 
partial decay widths are obtained by using Eq. with C33 = | then multiplying by 9'^. 
The expected widths of these states are of order 10 keV, but more precise estimates cannot 
be made until the masses of the states and the coupling g are known. For the nonstrange 
doubly heavy baryons, in the limit of infinite heavy quark mass, we obtain 

^.! = 55MeV r^Vr V 



TIE'] = TIE'*] = -^pi = 55 MeV {^] — ^ (34) 

^ ^ ^ ^ 47r/2-^" VO-5/ V250MeV/ ^ ^ 

16 



for the total widths, and for the branching fractions we find 

Br[S'* ^ H* ttI 5 BrE' ^ E* vrl 



(35) 



Br[S'*^S7r] 4 Br[S'^S7r] 
These relations receive large corrections due to phase space effects. Once the hyperfine 
splittings are taken into account the factors of p^. will differ greatly for the four decays. 
To get a feeling for these effects in the doubly charm sector we choose = 3440 MeV, 
m=. = 3520 MeV, m=/ = 3740 MeV, and m=*/ = 3820 MeV, which corresponds to a 
diquark excitation energy of 300 MeV and hyperfine splittings of 80 MeV. We then find 

r [E'J = f 336 MeV V [S*^] = f 78 MeV 

V\E'* ^ S* vrl V\E' —> S* vrl 

= 0.56 = 2.3 . (36) 

Note that the S'^^ unlike the S^* strongly prefers to decay to S*^ relative to Sec despite the 
phase space suppression. This may be useful for distinguishing S'^* and S'^^ experimentally. 

The SELEX S+"''(3780) is broad relative to the other SELEX doubly charm candidates. 
Since it is 260 MeV heavier than the S^(3520), it is a natural candidate for one of the low 
lying excited doubly charm baryons. Unfortunately, no measurement of the width exists 
and the pattern of decays is also hard to understand, since Ref. states that 50% of the 
decays to K^K^-k^-k^ are through S^(3520) vr^ while the other 50% are weak decays. More 
information on the quantum numbers of the S++(3780) and the H+(3520) are needed before 
we can determine which of the excited doubly charm baryons should be identified with the 
2++ (3780). 

In this paper we have developed a generalization of HH^PT which incorporates heavy 
quark-diquark symmetry and includes the leading symmetry breaking corrections from the 
chromomagnetic couplings of the heavy quark and diquark. We also included electromag- 
netic interactions in the Lagrangian, and obtained an estimate of the width of the J = \ 
member of the ground state doubly charm baryon doublet. The width of this state is com- 
pletely dominated by electromagnetic decays. We showed how to include the lowest lying 
doubly charm baryons which are expected to be excitations of the doubly charm diquark 
rather than the light degrees of freedom. Strong decay widths of low lying excited states were 
calculated and the states are expected to be rather narrow because of limited phase space 
available for the decays. Of particular interest is the doubly charm strange sector where we 
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expect three pairs of excited baryons whose strong decay must violate isopsin conservation 
because they are below the kaon decay threshold. These states will have narrow widths of 
10 keV or less. Experimental efforts to observe the narrow doubly charm strange baryons 
would be of great interest. 
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